Celestial Computing with Numerit 



Greenwich Apparent Sidereal Time 

This Numerit program (demogast) demonstrates how to calculate the Greenwich 
apparent sidereal time using both a low and high precision algorithm. Both numerical 
methods are based on the IAU 1980 nutation algorithm. 

function gastl 

This function calculates the apparent Greenwich sidereal time using the first few terms of 
the IAU 1980 nutation algorithm. 

The Greenwich apparent sidereal time is given by the expression 

9 = On, + Ay cos (e m + Ae) (1) 

In this equation 0 m is the Greenwich mean sidereal time, A\)/ is the nutation in longitude, 

£ m is the mean obliquity of the ecliptic and Ae is the nutation in obliquity. 

The Greenwich mean sidereal time, in degrees, is calculated using the expression 

9 m = 100.46061837 + 36000.770053608T + 0.000387933T 2 — TV38710000 (2) 

where T = (JD - 2451545)/36525 and JD is the Julian date on the UT1 time scale. 

The mean obliquity of the ecliptic is determined from 

£ m = 23°26 / 21. // 448 - 46."8150r- 0."00059T 2 + 0/TJ01813T 3 (3) 

The nutations in obliquity and longitude involve the following three trigonometric 
arguments (in degrees): 



L = 280.4665 + 36000.7698T 



L' = 218.3165 + 48 1267.88 13T 

Q. = 125.04452 - 1934. 13626 IT 
The calculation of the nutations use the following two equations: 

A\|/ = -17.20 sin Q. - 1.32 sin2L - 0.23 sin2Z/ + 0.21 sin 2 Q. 

Ae = 9.20cos£2 + 0.57cos2L + 0.10cos2Z/ - 0.09cos2f2 
where these corrections are in units of arc seconds. 



( 4 ) 
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The syntax of this function is as follows: 

function gastl (jdate, gst) 

Greenwich apparent sidereal time 
input 

jdate = Julian date 
output 

gst = Greenwich apparent sidereal time (radians) 

(0 <= gst <= 2 pi) 

function gast2 

This function calculates the mean or apparent Greenwich sidereal time. For the apparent 
sidereal time calculation, the obliquity in longitude and obliquity are determined using a 
Numerit function named nutation which implements the full IAU 1980 nutation 
algorithm (108 terms). Please note that the Julian date can be passed to this function in a 
high-order (integer) and low-order (fractional) part. In the argument list k determines the 
type of Greenwich sidereal time calculation (k = 0 => mean, k = 1 => apparent). This 
function was ported to Numerit using the Fortran version of the NOVAS (Naval 
Observatory Vector Astrometry Subroutines) source code which was developed at 
the United States Naval Observatory. 

For this numerical method, the nutation in longitude is determined from a trigonometric 
series of the form 



n 

A\)/ = ^ S j sin A i 

i= 1 



The nutation in obliquity is determined from a series of the form 



(6) 



n 

As = ^ C t cos A i 

i= 1 



(7) 



where 



A. = a j l + b j 1' + c t F + D + e t Cl 
In this expression /, f , F, D and Cl are fundamental arguments. 
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The syntax for this Numerit function is as follows: 

function gast2 (tjdh, tjdl, k, gst) 

this function computes the greenwich sidereal time 
(either mean or apparent) at julian date tjdh + tjdl 

input 

tjdh = julian date, high-order part 

tjdl = julian date, low-order part 

julian date may be split at any point, but for 
highest precision, set tjdh to be the integral part of 
the julian date, and set tjdl to be the fractional part 

k = time selection code 

set k=0 for greenwich mean sidereal time 
set k=l for greenwich apparent sidereal time 

output 

gst = greenwich (mean or apparent) sidereal time in hours 

The program will prompt you for the calendar date and universal time. The following is a 
typical draft output created with this program: 

program demogast - Greenwich apparent sidereal time 

calendar date October 21, 1999 

universal time 10 h 20 m 30 s 

Julian date 2451472.930902778 

low precision 12 h 18 m 11.32 s 

high precision 12 h 18 m 11.31 s 
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